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MULTIPATH DATA ACQUISITION 
SYSTEM AND METHOD 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This application is related to U.S. Application Serial No. , filed on even 

5 date herewith, by Randy K. Roushall and Robert K. Crawford, and entitled 
"Phase-Shifted Data Acquisition System and Method," which is incorporated 
herein by reference. 

TECHNICAL FIELD 

This invention relates to data acquisition systems and methods. 

10 BACKGROUND 

Data acquisition systems and methods may be used in a variety of 
applications. For example, data acquisition techniques may be used in nuclear 
magnetic resonance imaging systems and Fourier transform spectrometer systems. 
Such techniques also may be used in mass spectrometer systems, which may be 

15 configured to determine the concentrations of various molecules in a sample. A 
mass spectrometer operates by ionizing electrically neutral molecules in the 
sample and directing the ionized molecules toward an ion detector. In response 
to applied electric and magnetic fields, the ionized molecules become spatially 
separated along the flight path to the ion detector in accordance with their mass- 

20 to-charge ratios. 

Mass spectrometers may employ a variety of techniques to distinguish ions 
based on their mass-to-charge ratios. For example, magnetic sector mass 
spectrometers separate ions of equal energy based on their momentum changes in 
a magnetic field. Quadrupole mass spectrometers separate ions based on their 

25 paths in a high frequency electromagnetic field. Ion cyclotrons and ion trap mass 
spectrometers distinguish ions based on the frequencies of their resonant motions 
or stabilities of their paths in alternating voltage fields. Time-of-flight (or "TOF") 
mass spectrometers discriminate ions based on the velocities of ions of equal 
energy as they travel over a fixed distance to a detector. 

30 In a time-of-flight mass spectrometer, neutral molecules of a sample are 

ionized, and a packet (or bundle) of ions is synchronously extracted with a short 
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voltage pulse. The ions within the ion source extraction are accelerated to a 
constant energy and then are directed along a field-free region of the 
spectrometer. As the ions drift down the field-free region, they separate from one 
another based on their respective velocities. In response to each ion packet 
5 received, the detector produces a data signal (or transient) from which the 

quantities and mass-to-charge ratios of ions contained in the ion packet may be 
determined. In particular, the times of flight between extraction and detection 
may be used to determine the mass-to-charge ratios of the detected ions, and the 
magnitudes of the peaks in each transient may be used to determine the number 

10 of ions of each mass-to-charge in the transient. 

A data acquisition system (e.g., an integrating transient recorder) may be 
used to capture information about each ion source extraction. In one such 
system, successive transients are sampled and the samples are summed to 
produce a summation, which may be transformed directly into an ion intensity 

15 versus mass-to-charge ratio plot, which is commonly referred to as a spectrum. 
Typically, ion packets travel through a time-of-flight spectrometer in a short time 
(e.g., 100 microseconds) and ten thousand or more spectra may be summed to 
achieve a spectrum with a desired signal-to-noise ratio and a desired dynamic 
range. Consequently, desirable time-of-flight mass spectrometer systems include 

20 data acquisition systems that operate at a high processing frequency and have a 
high dynamic range. 

In one data acquisition method, which has been used in high-speed digital- 
to-analog converters, data is accumulated in two or more parallel processing 
channels (or paths) to achieve a high processing frequency (e.g., greater than 100 

25 MHz). In accordance with this method, successive samples of a waveform (or 
transient) are directed sequentially to each of a set of two or more processing 
channels. The operating frequency of the components of each processing channel 
may be reduced from the sampling frequency by a factor of N, where N is the 
number of processing channels. The processing results may be stored or 

30 combined into a sequential data stream at the original sampling rate. 
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SUMMARY 

When applied to applications in which sample sets (or transients) are 
accumulated to build up a composite signal (e.g., TOF mass spectrometer 
applications), the process of accumulating samples in parallel processing channels 
5 may introduce noise artifacts that are not reduced by summing the samples from 
each processing channel. In particular, although contributions from random noise 
and shot noise may be reduced by increasing the number of transients summed, 
each processing channel may contribute to the composite signal a non-random 
pattern noise that increases with the number of transients summed. Such pattern 

10 noise may result from minute differences in digital noise signatures induced in the 
system by the different parallel processing paths. For example, the physical 
separations between the components (e.g., discrete memory, adders and control 
logic) of a multi-path or parallel-channel data acquisition system may generate 
voltage and current transitions within the board or chip on which the data 

15 acquisition system is implemented. The unique arrangement of each processing 
path may induce a unique digital noise signature (or pattern noise) in the analog 
portion of the system. The resulting digital noise signature increases as the 
composite signal is accumulated, limiting the ability to resolve low-level transient 
signals in the composite signal. 

20 The invention features improved data acquisition systems and methods 

that substantially reduce accumulated pattern noise to enable large numbers of 
data samples to be accumulated rapidly with low noise and high resolution. 

In one aspect of the invention, a data acquisition system includes an 
accumulator that has two or more parallel accumulation paths and accumulates 

25 corresponding data samples across a transient sequence through different 
accumulation paths. 

As used herein, the phrase "corresponding data samples across a transient 
sequence" refers to the summation of data samples from different transients 
having similar mass-to-charge ratios. 

30 Embodiments may include one or more of the following features. 

A controller preferably is coupled to the accumulator and preferably is 
configured to cycle the accumulation of data samples through each of the 
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accumulation paths. The controller preferably is configured to selectively enable 
each accumulation path. 

Each accumulation path may include an adder and a memory. The 
accumulation path memory may comprise a dual port random access memory. 
5 Each accumulation path preferably is configured to produce an output 

representative of the sum of two inputs. The accumulation paths may be coupled 
in series with a first input of each accumulation path coupled to the sampler and 
a second input of each accumulation path coupled to the output of another 
accumulation path. 
10 The data acquisition system may include an ion detector. 

In another aspect, the invention features a time-of-flight mass 
spectrometer that includes an ion detector, a sampler, and an accumulator. The 
ion detector is configured to produce a transient sequence from a plurality of ion 
packets. The sampler is configured to produce a plurality of data samples from 
15 the transient sequence. The accumulator comprises two or more accumulation 
paths and accumulates corresponding data samples across the transient sequence 
through different accumulation paths. 

In another aspect, the invention features a method of acquiring data. In 
accordance with this inventive method, a plurality of data samples is produced 
20 from a transient sequence, and corresponding data samples are accumulated 
across the transient sequence through two or more parallel accumulation paths. 
Embodiments may include one or more of the following features. 
The accumulation of data samples preferably is cycled through each of the 
parallel accumulation paths. The data samples may be cycled by selectively 
25 enabling each accumulation path. Alternatively, the data samples may be cycled 
by selectively directing consecutive data sample sets to a respective accumulation 
path. An analog transient may be converted into one or more digital data 
samples. A transient may be produced from a received ion packet. A plurality of 
packets may be launched along a flight path defined in a time-of-flight mass 
30 spectrometer. 

Among the advantages of the invention are the following. 
By accumulating corresponding data samples across a transient sequence 
through different accumulation paths, the overall noise level induced in the 
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spectrum data may be reduced. This feature improves the signal-to-noise ratio in 
the resulting spectrum and, ultimately, improves the sensitivity of the data 
acquisition system. 

Other features and advantages of the invention will become apparent from 
5 the following description, including the drawings and the claims. 

DESCRIPTION OF DRAWINGS 

FIG. 1 is a block diagram of a time-of-flight mass spectrometer, including a 
flight tube and a data acquisition system. 

FIG. 2A is a plot of a transient sequence produced by an ion detector in the 
10 flight tube of FIG. 1. 

FIG. 2B is a diagrammatic view of a plurality of sets of data samples 
produced by the data acquisition system from transient sequence of FIG. 2 A. 

FIG. 2C is a plot of an accumulated sample spectrum produced by the data 
acquisition system from the data sample sets of FIG. 2B. 
15 FIG. 2D is a diagrammatic view of an accumulated data sample set 

corresponding to the accumulated sample spectrum of FIG. 2C. 

FIG. 3 is a block diagram of the data acquisition system of FIG. 1, including 
a plurality of accumulation paths each having a respective accumulator. 

FIG. 4 is a block diagram of an accumulator of the data acquisition system 
20 of FIG. 3. 

FIG. 5 is a plot of signals of a mass spectrometer having a single-path 
accumulator that is clocked by an accumulation clock that is synchronized with a 
sampling clock. 

FIG. 6 is a plot of signals of a mass spectrometer having a single-path 
25 accumulator that is clocked by an accumulation clock that is shifted in phase 
relative to a sampling clock. 

FIG. 7 is a plot of signals of a mass spectrometer having a accumulator 
with multiple accumulation paths, each of which is clocked by a respective 
accumulation clock that is shifted in phase relative to a sampling clock by a 
30 respective amount. 
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DETAILED DESCRIPTION 

Referring to FIG. 1, a time-of- flight mass spectrometer 10 includes an ion 
source 12, a flight tube 16, a data acquisition system 18, and a processor 20 (e.g., 
a computer system). Time-of-flight mass spectrometer 10 may be arranged in an 
5 orthogonal configuration or on-axis configuration. Ion source 12 may generate 
ions using any one of a variety of mechanisms, including electron impact, 
chemical ionization, atmospheric pressure ionization, glow discharge and plasma 
processes. Flight tube 16 includes an ion detector 22 (e.g., an electron 
multiplier), which is configured to produce a sequence of transients 24 containing 

10 a series of pulses from which the quantities and mass-to-charge ratios of the ions 
within each transient may be determined. In operation, sample molecules are 
introduced into source 12, ion source 12 ionizes the sample molecules, and 
packets of ionized molecules are launched down flight tube 16. A conventional 
orthogonal pulsing technique may be used to release the packets of ions into flight 

15 tube 16. The ions of each packet drift along a field-free region defined inside 
flight tube 16. As they drift down flight tube 16, the ions separate spatially in 
accordance with their respective masses, with the lighter ions acquiring higher 
velocities than the heavier ions. In FIG. 1, an ion packet 26 consists of two 
constituent ion concentrations: a relatively low concentration of lighter ions 28, 

20 and a relatively high concentration of heavier ions 30. 

Referring to FIGS. 2A-2D, after an initial time delay corresponding to the 
time between the extraction pulse and the arrival of the first (i.e., the lightest) 
ions at the detector, detector 22 produces a transient 32 representative of the ion 
intensities in the detected ion source extraction. The peaks 34, 36 of transient 32 

25 represent the numbers of light ions 28 and heavy ions 30, respectively, and the 
peak times correspond to the mass-to-charge ratios of the ions within transient 32. 
Detector 22 produces a sequence of additional transients 38, 40 from subsequent 
ion packets launched into flight tube 16. Data acquisition system 18 samples m 
transients 32, 38, 40, and produces from each transient data samples (d jtl , d. 2 , . . . 

30 > d jim , where j = 1 to k) that may be represented as a respective data sample set 
42, 44, 46 (FIG. 2B). The resulting data samples (d M , d. 2 , . . . , d f J are 
accumulated by data acquisition system 18 to produce a spectrum 48 (FIG. 2C), 
which may be represented by an accumulated data sample set 50 (FIG. 2D), in 
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which each member corresponds to the sum of ion samples (d hi , where i = 1 to 
m) having similar mass-to-charge ratios. 

Data acquisition system 18 may be designed to control the operation of 
time-of- flight mass spectrometer 10, collect and process data signals received from 
5 detector 22, control the gain settings of the output of ion detector 22, and provide 
a set of time array data to processor 20. As explained in detail below, data 
acquisition system 18 is configured to accumulate corresponding data samples 
across the transient sequence 24 through each of a plurality of parallel data 
accumulation paths. In this way, data acquisition system 18 may accumulate 

10 data samples at a high speed, while reducing the impact of noise introduced by 
data acquisition system 18. 

Referring to FIG. 3, in one embodiment, data acquisition system 18 
includes a sampler 60 (e.g., a high speed flash analog-to-digital converter), a 
multipath sample accumulator 62 and a controller 64. Sampler 60 samples 

15 transients 24 and produces a series of data samples 65, which are applied to an 
input of sample accumulator 62. The output of sampler 60 is a series of digital 
signals (i.e., an n-bit word) each of which represents instantaneous ion intensities 
at respective sampling times. The resolution with which sampler 60 captures the 
instantaneous ion intensities is determined by the bit width of sampler 60. 

20 Sample accumulator 62 includes a plurality (N) of accumulators 66 that define a 
respective plurality of parallel data accumulation paths. In operation, controller 
64 directs the data samples to one of the N accumulators 66 in sequence. Thus, 
each accumulator 66 processes only 1/N of the data samples and need only 
operate at a frequency that is roughly only 1/N of the operating frequency of a 

25 comparable single-path data acquisition system (e.g., the sampling rate). At the 
same time, controller 64 cycles the accumulation of data samples through each of 
the accumulation paths so that corresponding data samples across the transient 
sequence are accumulated through each of the accumulation paths. For example, 
assuming that eight data samples (d M> d 2 . . . , d 8 ) are measured for each 

30 transient 24, the data samples would be accumulated after each of m transients as 
follows: 
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Table 1 : Cycled Transient Accumulation 



As explained in detail below, each accumulation path induces a unique 
5 noise signal in each of the transients 24. By cycling the accumulation of data 
samples through each of the N accumulation paths, data acquisition system 18 
reduces the noise level in the accumulated spectrum 48 relative to a system that 
does not perform such cycling. In particular, the accumulated spectrum may be 
expressed as: 

10 

D(h) = d ( h > i) CD 

where d(h, j) is the j lh accumulated data point having a mass-to-charge ratio of h. 
The component data samples of the accumulated data points (d(h, j)) may be 
15 expressed as follows: 
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d(h, j) = s(h, j) + v(h, j) + n(h, j) 



(2) 



where s(h, j) is the noise-free signal, v(h, j) is the signature (or pattern) noise 
induced by the paths of the data acquisition system, and n(h, j) is random noise. 

5 The induced signature noise (v(h, j)) is a non-random, non-white noise source 
that is specific to each accumulation path. In a dual-path data accumulation 
embodiment, all of the even-numbered samples have the same induced digital 
noise (i.e., v(2, j) = v(4, j)), and all of the odd-numbered samples have the same 
induced digital noise (i.e., v(l, j) = v(3, j)). Similarly, for a four-path data 

10 accumulation embodiment, v(l, j) = v(5, j), v(2, j) = v(6, j), v(3, j) = v(7, 
j), andv(4, j) = v(8, j). 

Without path cycling, the induced signature noise is the same across the 
data samples (i.e., v(h, 1) = v(h, 2) = . . . = v(h, m)). As a result, the 
accumulated spectrum signal may be estimated by the following equation: 



The random noise source (n(h, j)) falls off by the square root of m and, therefore, 
becomes negligible for large values of m. The induced signature noise (v(h)), 
20 however, increases because it is specific to each an accumulation channel and not 
random. Thus, in a dual-path data accumulation system, 



For large transient* signals, the s(h) term dominates the v(h) and, consequently, 
the data acquisition system may resolve the data signal. For small transient 
signals, however, the v(h) term may be larger than the s(h) term, making it 
difficult to resolve the data signal. In particular, for small transient signals, the 
30 difference between data points in the accumulated spectrum may be estimated as 
follows: 



15 



D(h) = m * s(h) + m • v(h) + 2% n(h, j) 



(3) 



D(l) = m • s(l) + m • v(l) 
D(2) = m • s(2) + m * v(2) 



(4) 
(5) 



D(2) - D(l) = m • v(2) - m * v(l) 



(6) 
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This difference is the cause of the induced pattern noise signal 94 shown in FIG. 
6. 

On the other hand, if the sample accumulation is cycled through each of 
5 the N accumulation paths as described above, the induced digital noise signatures 
may be reduced substantially or eliminated as follows. In a dual-path data 
accumulation embodiment the following relationships are established (ignoring 
random noise). The data samples for the first transient may be expressed as 
follows: 



d(l, 1) 
d(2, 1) 
d(3, 1) 
d(4, 1) 



d(l, 2) 

20 d(2, 2) 

d(3, 2) 
d(4, 2) 



= s(l, 1) + v(l, 1) (7) 

= s(2, 1) + v(2, 1) (8) 

= s(3, 1) + v(l, 1) (9) 

= s(4, 1) + v(2, 1) (10) 



= s(l, 2) + v(2, 2) (11) 

= s(2, 2) + v(l, 2) (12) 

= s(3, 2) + v(2, 2) (13) 

= s(4, 2) + v(l, 2) (14) 



where v(l, 1) = v(3, 1) and v(2, 1) = v(4, 1) in a dual-path data accumulation 
system. The data samples for the second transient may be expressed as follows: 



Since the induced digital signature noise (v(h, j) is the same for all transients (i.e., 
25 v(l, 1) = v(l, 2) and v(2, 1) = v(2, 2)), equations (11)-(14) may be re-written as 
follows: 



30 



d(l, 2) 
d(2, 2) 
d(3, 2) 
d(4, 2) 



= s 



= s 



(1, 2) + v(2, 1) 

s(2, 2) + v(l, 1) 

s(3, 2) + v(2, 1) 

;(4, 2) + v(l, 1) 



(15) 
(16) 
(17) 
(18) 
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Thus, the summation of the data points for the first two transients may be 
expressed as follows: 



D(l) = s(l, 1) + s(l, 2) + [v(l, 1) + v(2, 1)] (19) 

5 D(2) = s(2, 1) + s(2, 2) + [v(2, 1) + v(l, 1)] (20) 

D(3) = s(3, 1) + s(3, 2) + [v(l, 1) + v(2, 1)] (21) 

D(4) = s(4, 1) + s(4, 2) + [v(2, 1) + v(l, 1)] (22) 



As a result, the induced digital signature noise terms drop out in the difference 
10 between any two adjacent data points. For example, the difference between the 
first accumulated data point (D(l)) and the second accumulated data point (D(2)) 
may be expressed as follows: 

D(2) - D(l) = [s(2, 1) + s(2, 2)] - [s(l, 1) + s(l, 2)] (23) 

15 

In general, the difference between any two adjacent data points may be expressed 
as follows: 

D(h) - D(h-l) = ^ [s(h, j) + s(h-l, j)] + L",., [n(h, j) + n(h-l, j)] (24) 

20 

The only noise term remaining in equation (24) is the random noise source (n(h, 
j)), which drops off by the square root of the number of summations (m). In this 
case, equation (3) reduces to the following form: 

25 D(h) = m • s(h) + X 11 ^ n(h, j) (25) 

This feature of the data acquisition system advantageously improves the signal-to- 
noise ratio of the accumulated spectrum 48 and, ultimately, improves the 
30 sensitivity of the measurements of mass spectrometer 10. 

Referring to FIG. 4, in one embodiment, each accumulator 66 includes an 
adder 70 and a memory system 72. In operation, during each clock cycle adder 
70 computes the sum of the signal values applied to inputs 74, 76, and memory 
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system 72 stores the computed sum. As shown in FIG. 4, memory system 72 may 
include an input address counter 78, an output address counter 80 and a dual port 
random access memory (RAM) 82. In one embodiment, controller 64 selectively 
enables adder 70 so that corresponding data samples generated by sampler 60 are 
5 accumulated through each of the data accumulation paths. In another 
embodiment, controller 64 selectively directs data samples to respective 
accumulation paths, for example, by controlling the output of a 1-by-N 
multiplexer, which is coupled between sampler 60 and sample accumulator 62. 
Other embodiments are within the scope of the claims. 

10 Referring to FIG. 5, in a single accumulation path embodiment, sampler 60 

is configured to sample transients 24 received from ion detector 22 in response to 
the falling edge of a sampling clock 90. Sample accumulator 62, on the other 
hand, is configured to accumulate data in response to the rising edge of an 
accumulation clock 92. If sampling clock 90 and accumulation clock 92 are in 

15 phase (as shown), the rising edge of accumulation clock 92 may induce a noise 
signal 94 in an analog transient 98. The induced noise ultimately may appear in 
data samples 96 produced by sampler 60, reducing the signal-to-noise ratio and 
reducing the sensitivity of the accumulated spectrum 48. Without being limited to 
a particular theory, it is believed that this noise is generated, at least in part, by a 

20 capacitive coupling between sample accumulator 62 and sampler 60. 

The magnitude of the accumulation clock induced noise signal 94 may be 
reduced substantially by shifting the phase of accumulation clock 92 relative to 
sampling clock 90. For example, referring to FIG. 6, by shifting accumulation 
clock 92 relative to sampling clock 90, the noise signal peaks 99, which are 

25 induced in transient 98, may be shifted away from the sampling times (i.e., the 
falling edges of sampling clock 90) to reduce the noise level appearing in 
accumulated spectrum 48. Accumulation clock 92 preferably is shifted relative to 
sampling clock 90 by an amount selected to minimize induced noise signal 94. In 
one embodiment, accumulation clock 92 preferably is shifted between 90° and 

30 270° relative to sampling clock 90, and more preferably is shifted approximately 
180° relative to sampling clock 90. 

Referring to FIG. 7, in another embodiment, sample accumulator 62 
includes two accumulation paths (Path A and Path B), each of which accumulates 



12 



data samples in response to a respective accumulation clock 100, 102. In this 
embodiment, the phase of each accumulation clock 100, 102 is shifted relative to 
sampling clock 90 by a respective amount selected to reduce the overall noise in 
the accumulated spectrum 48. The phases of accumulation clocks 100, 102 may 
5 be shifted by the same amount relative to sampling clock 90, or they may be 
shifted independently by different amounts (as shown) . 

The above-described phase shift between sampling clock 90 and the one or 
more accumulation clocks may be implemented by a multiphase frequency 
synthesizer 110 (FIG. 3) that includes a phase-locked loop, a delay-locked loop, or 

10 any phase-shifting clock driver. In addition, the phase shift between sampling 
clock 90 and the one or more accumulation clocks may be programmable to 
enable the relative clock phases to be adjusted during an initial calibration of 
mass spectrometer 10 or dynamically during operation of mass spectrometer 10. 

The systems and methods described herein are not limited to any particular 

15 hardware or software configuration, but rather they may be implemented in any 
computing or processing environment. Data acquisition controller 64 preferably 
is implemented in hardware or firmware. Alternatively, controller 64 may be 
implemented in a high level procedural or object oriented programming language, 
or in assembly or machine language; in any case, the programming language may 

20 be a compiled or interpreted language. 

Still other embodiments are within the scope of the claims. 
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